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Abstract We review the influence of a magnetic background on the phase diagram 
of strong interactions and how the chiral and deconfining transitions can be affected. 
First we summarize results for both transitions obtained in the framework of the 
^ C| linear sigma model coupled to quarks and to the Polyakov loop, and how they com- 

Q_[ pare to other effective model approaches and to lattice QCD. Then we discuss the 

outcome of the magnetic MIT bag model that yields a behavior for the critical de- 
_ confining temperature which is compatible with recent lattice results and magnetic 

_C catalysis. The qualitative success of the magnetic MIT bag model hints to T c being 

a confinement-driven quantity, and leads us to the discussion of its behavior as pre- 
dicted within the large-A^ t limit of QCD, which is also in line with the most recent 
lattice QCD results provided that quarks behave paramagnetically. 



1 Introduction 

The thermodynamics of strong interactions under a strong magnetic background has 
proven to be a very rich and subtle subject. Recent developments were initially mo- 
tivated by the utility of magnetic fields in separating charge in space, which would 
render the possible formation of sphaleron-induced CP-odd domains in the plasma 
created in high-energy heavy ion collisions, in the so-called chiral magnetic ef- 
fect HI, measurable. In fact, the magnetic fields created in non-central collisions 
in heavy ion experiments at RHIC-BNL and the LHC-CERN are possibly the high- 
est since the epoch of the electroweak phase transition, reaching values such as 
B ~ 10 19 Gauss (eB ~ 6m^) for peripheral collisions at RHIC 12) and even much 
higher at the LHC due to the fluctuations in the distribution of protons inside the 
nuclei [3|. 
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From the theoretical point of view, the non-trivial role played by magnetic fields 
in the nature of phase transitions has been known for a long time (4). Modifica- 
tions in the vacuum of QED and QCD have also been investigated within different 
frameworks, mainly using effective models |!5, S] [7j [8] |9] [lOl [XI] [T2) , especially the 
NJL model ff3l . and chiral perturbation theory (14\ [T5l [T6l . but also resorting to 
the quark model ifTTl and certain limits of QCD fl8l . Interesting phases in dense 
systems 11911201121112211 . as well as effects on the dynamical quark mass H231 were 
also considered. Nevertheless, the mapping of the new T — eB phase diagram is still 
an open problem. There are clear indications that sufficiently large magnetic fields 
could significantly modify the behavior of the chiral and the deconfinement phase 
transition lines ||24l |231l26ll271l2^|29ll30l[3T1l32l[33l |34l |35ll36ll37ll38ll39ll40ll4Tl, 
or even transform the vacuum into a superconducting medium via p -meson conden- 
sation [42 1 . Although most of the analyses so far relied on effective models, lattice 
QCD has definitely entered the field and has been producing its first results for the 
phase diagram |43ll44ll45ll46ll47l. 
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Fig. 1 Originally expected magnetic field-temperature phase diagram of strong interactions. The 
thick lines indicate first-order transitions, the filled circles are the (second-order) endpoints of these 
lines, and the thin dashed lines stand for the corresponding crossovers. A new phase with broken 
chiral symmetry and deconfinement appears at high magnetic fields. Extracted from Ref. 1271 . 



From the first results obtained within effective models for the deconfining l24l 
and chiral (25 1 transition lines, one would expect the phase diagram structure illus- 
trated in Fig. [T] as discussed in Ref. 1271 . Indeed, after the prediction of a split- 
ting between the chiral and deconfining transition lines, with the appearance of 
a new phase, in Ref. l27l . several model descriptions produced the same effect 
EOl [3T] [32] [33 [39]|. However, until 2011, all model studies have yielded either a 
monotonically increasing or an essentially flat functional form for the deconfine- 
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ment critical line as B increases to very large valuefl Pioneering lattice simulations 
l44l also found an essentially flat behavior for both critical lines, that seemed to 
increase together at a very low rate. Nevertheless, since the pion mass used in these 
simulations was still very high, this could be an indication that one would proba- 
bly need huge magnetic fields in the simulations in order to be able to compare to 
effective model predictions. 

This was the scenario, rather coherent in terms of expectations for the behavior of 
the critical lines for the chiral and deconfining transitions in the presence of a mag- 
netic background, until lattice simulations of magnetic QCD with physical masses 
and fine grids were performed [46} and showed that both critical temperatures actu- 
ally go down for increasing B, saturating for very large fields, very differently from 
what has been predicted by all previous effective model calculations and found in 
previous lattice simulations. 

Soon after the appearance of the new lattice results, the behavior of the critical 
temperature for deconfinement in the presence of a very large magnetic field was 
addressed within the MIT bag model BUI , a very economic model in terms of pa- 
rameters to be fixed (essentially one) and other ingredients usually hard to control 
in more sophisticated effective theories. The model is, of course, crude in numerical 
precision and misses the correct nature of the (crossover) transition. Nevertheless, 
it provides a simple setup for the discussion of some subtleties of vacuum and ther- 
mal contributions in each phase. It was shown in Ref. |40]| that the influence of the 
magnetic field on the thermodynamics of both extreme energy domains is captured, 
so that the model furnishes a reasonable qualitative description of the behavior of 
the critical temperature in the presence of B, decreasing and saturating. 

The fact that chiral models, even when coupled with the (static) Polyakov loop 
sector, seem to fail in the description of the behavior of T c x eB, whereas the (as- 
sumedly simple) MIT bag approach finds a good qualitative agreement, suggests 
that the critical temperature in QCD is a confinement-driven observable. This was 
also hinted by a previous successful description of the behavior of the critical tem- 
perature as a function of the pion mass and isospin chemical potential, as compared 
to lattice data, where chiral models also failed even qualitatively [48 49] . If con- 
finement dynamics plays a central role in guiding the functional behavior of T c , a 
the large-A^ limit of QCD should provide an adequate and powerful framework to 
study associated magnetic thermodynamics. In fact, it was shown in Ref. ||5Q| on 
very general grounds that the fact that the deconfining temperature decreases and 
tends to saturate for large B, although this last point cannot be proven in a model- 
independent way, depends solely on quarks behaving paramagnetic ally. 

In the sequel we summarize results for the chiral and deconfining transitions ob- 
tained in the framework of the linear sigma model coupled to quarks and to the 
Polyakov loop, especially the prediction of a splitting of the two critical lines, and 
how they compare to other effective model approaches as well as to lattice QCD. 



1 Contrastingly, a significant decrease in the critical temperature as a function of B, vanishing at 
eB c ~ 25m 2 K , was found in Ref. 1241 . featuring the disappearance of the confined phase at large 
magnetic fields. This phenomenon that was not reproduced by any other effective model nor ob- 
served on the lattice (even for much larger fields). 
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Then we discuss the outcome of the magnetic MIT bag model that yields a behav- 
ior for the critical deconfining temperature compatible with the most recent lattice 
simulations and magnetic catalysis. We continue with a discussion of very recent re- 
sults, starting with the rather general findings within the large-A^ c limit of magnetic 
QCD. Finally, we present our conclusions. 



2 Modified dispersion relations and integral measures 

In the presence of a classical, constant and uniform (Abelian) magnetic field, dis- 
persion relations and momentum integrals will be modified. In order to compute 
vacuum and thermal determinants and Feynman diagrams, it is necessary to express 
these quantities in a convenient fashion. Lorentz invariance is broken by the pre- 
ferred direction established by the external field, and Landau orbits redefine the new 
counting of quantum states J4) . 

For definiteness, let us take the direction of the magnetic field as the z-direction, 
B = Bz. One can compute, for instance, the modified effective potential or the modi- 
fied pressure to lowest order by redefining the dispersion relations of charged scalar 
and spinor fields in the presence of B, using the minimal coupling shift in the gradi- 
ent and the field equations of motion^. For this purpose, it is convenient to choose 
the gauge such that = (A , A) = (0, -By, 0,0). 

For scalar fields with electric charge q, such as pions, one has 

{d 2 + m 2 )<p=0, (1) 
d^^d^+iqA^. (2) 

After decomposing into Fourier modes, except for the dependence in the coordi- 
nate y, one obtains 



<p"{y) + 2m 



p\-p\-rn 2 \ q 2 B 2 ( pA 2 ' 



2m ) 2m \^ ^ qB J 



<P0)=0, (3) 



which has the form of a Schrodinger equation for a harmonic oscillator. Its eigen- 
modes correspond to the well-known Landau levels 

where <x>b = \q\Bjm and £ is a positive ((, > 0) integer, and provide the new disper- 
sion relation: 

p 2 0n =p 2 + m 2 + (2£+l)\q\B. (5) 



2 Higher-order (loop) corrections need the full propagator, not only its poles. 
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One can proceed in an analogous way for fermions with charge q. From the free 
Dirac equation {iy^d^ — m) y = 0, and the shift in <9 M , one arrives at the following 
Schrodinger equation 



u s (y) + 2m 



2 2 2 

Pa - Pz - m 



2m 



2r>2 



q\Bs\ q z B 



2m 



y+ 



qBj 



which yields the new dispersion relation for quarks: 

pl n =p 2 z +m 2 + (2£ + l-s)\q\B, 



u,(y)=0, (6) 



(7) 



where s = ±1 is the spin projection in the z direction. 

It is also straightforward to show that integrals over four momenta and thermal 
sum-integrals acquire the following forms, respectively |25ll5ll : 



/ 



d 4 k \q\B " 
2% '-' 



(2n)* 
d 3 k 



dko dk z 
2% 2k 

dk- 



2k 



n (=0 



2K 



(8) 
(9) 



where I represents the different Landau levels and n stands for the Matsubara fre- 
quency indices 



3 PLSM g effective model and the splitting of the chiral and 
deconfining transition lines 

Let us consider the two-flavor linear sigma model coupled to quarks and to the 
Polyakov loop, the PLSM (/ effective model, in the presence of an external magnetic 
field <27l . 

The confining properties of QCD are encoded in the complex-valued Polyakov 
loop variable L. As a matter of fact, the Polyakov loop sector only provides a de- 
scription of the behavior of the approximate order parameter for the Z(3) symmetry, 
which is explicitly broken by the presence of quarks. It is convenient for model- 
ing the deconfining transition and has a good agreement with lattice results for most 
thermodynamic quantities such as the pressure and energy density, especially for the 
pure glue theory, but it does not provide a dynamical description of confinement. 

The expectation value of the Polyakov loop L is an exact order parameter for 
color confinement in the limit of infinitely massive quarks: 



3 This will be a key feature in the discussion of recent results for the critical temperature, since T c 
seems to be a confinement-driven observable for both QCD transitions. 
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1/T 

Confinement: | j£j = ° ^ , L(x) = ^Tr^expJ/ 1 drA 4 (x,T)] , (10) 

o 

where Aa, — iAq is the matrix-valued temporal component of the Euclidean gauge 
field Af! and the symbol 3? denotes path ordering. The integration takes place over 
compactified imaginary time T, with periodic boundary conditions. 

The chiral features of the model are encoded in the dynamics of the 0(4) chiral 
field, which is an exact order parameter in the chiral limit, in which quarks and pions 
are massless degrees of freedom: 

Chiral svmmetrv . f (<*) ^ , low T (j> = (a, x) , 

Chiral symmetry • j (ty) = ; high T , n= (7r + ^ n -y ( n ) 

Here n is the isotriplet of the pseudoscalar pion fields and a is the chiral scalar field 
which plays the role of an approximate order parameter of the chiral transition in 
QCD, since chiral symmetry is explicitly broken by the nonzero quark masses. 

Within this effective model, the quark field y/ connects the Polyakov loop L and 
the chiral field (j>, making a bridge between confining and chiral properties. Quarks 
are also coupled to the external magnetic field since the u and d quarks are electri- 
cally charged. Thus, it is clear that the external magnetic field will affect the chiral 
dynamics as well as the confining properties of the model, as much as the latter can 
be captured by the Polyakov loop sector. 

This represents a natural generalization of the linear sigma model coupled to 
quarks l53l . an effective theory that has been widely used to describe different as- 
pects of the chiral transition, such as thermodynamic properties ll54l l55l l56l l57l l58l 
EH |60] |6T| IMl |63l IMl and the nonequilibrium phase conversion process 16511 . This 
generalization differs from previous ones l66l l67l by the inclusion of a bridge via 
the covariant derivative and, of course, because of the modifications brought about 
by the magnetic field. 

The Lagrangian of PLSM 9 describes the constituent quarks y/, which interact 
with the meson fields d, = (n 1 ± in 2 ) / \/2 and 7T° = 7T 3 , the Abelian gauge field 
a pi = ( fl °i a ) = (0, —By, 0,0), and the SI/ (3) gauge field A^ via the covariant deriva- 
tive = (d^ — iQa^ — ;A|,) with the charge matrix Q = diag(+2e/3, — e/3). Its 
explicit form is given by 

>: 

+ \D { *> | 2 - V (a, tt) - V L (L, T) , (12) 

(it) 

where = + iea^ is the covariant derivative acting on colorless pions. 
The chiral potential has the form 

V^o,n) = j(o 2 + n 2 -v 2 ) 2 -ho, (13) 



Jgf = 7 ]i?Df -g(a + iy 5 T ■7z)\ W + 1 - [(d^a) 2 + (d^ ) 2 ] 



Thermal chiral and deconfining transitions in the presence of a magnetic background 



7 



where h = f n m\, v 2 = fl - m % JX, A = 20, /„ ps 93MeV and m K ps 138MeV. The 
constituent quark mass is given by m q = m q ({o)) = g(o), and, choosing g = 3.3 
at T = 0, one obtains for the constituent quarks in the vacuum m q ps 310 MeV. At 
low temperatures quarks are not excited, and the model reproduces results from the 
usual linear <7-model without quarks. 

The Polyakov potential adopted is given by 

V L (L,T) L*L 2 ^r o y 



T 4 2 i=o 

V 3 



l-6L*L + A[L*" +U) - 3(L*L) 



(14) 



where Tq = T su ^ = 270MeV is the critical temperature in the pure gauge case and 
ao = 16 7T 2 /45 ps 3.51, a\ = -2.47, a 2 = 15.2, and b 3 = -1.75. Below we follow a 
mean-field analysis in which the mesonic sector is treated classically whereas quarks 
represent fast degrees of freedom. 





Fig. 2 Effects of temperature and magnetic field on quark confinenemt: The Polyakov loop poten- 
tial at T = 0.8To (top) and T = 1.27j) (bottom) and at zero magnetic field (left) and at eB = 9T 2 
(right). Extracted from Ref. 1271 . 



The one-loop corrections to the free energy Q coming from quarks can be written 



as: 



Jv 3d n q /T 



det(iy»D$ - m q ) 
dst^yv-dn-tnq) 



det r (>y D$ -m q 
dstlifDf -m q ) 



(15) 



so that the expectation values of the condensates can be obtained by minimizing the 
free energy 
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Q(a,L;T,B) = V$(o,n) +V L {L,T) + Q q (a,L,T) , 



(16) 



at fixed values of temperature and magnetic field. The interaction piece Q q (a,L 1 T) 
can be split into a vacuum (temperature-independent but still magnetic-field depen- 
dent) contribution and a thermal correction. The vacuum term has the form 



N c 



a%°(B) = —± £ \q f \B 



f=u,d 



£4 1) K) 

.n=t 



4 1; K) 



minus the standard vacuum correction in the absence of the magnetic field, 



2#c £ * B 

f=u,d 



(17) 



where we have defined the integral 



d d p 
J2k) u 



Vp 2 +m 2 



(18) 



and Mjj- = nty + 2l\qj\B. The thermal (paramagnetic) contribution is given by ||27l 

k 



para _ 



s= ±ie=ok=i 



where fl s £ is the energy of the i ,h Landau level at zero longitudinal momentum, 



Tr# A 



(19) 



1 /2 

(2£+ 1 — 2i)|^|Z?] , and the untraced Poly akov loop is such that 



Re [Tr^^] = Y?i=\ cos(k(pi), the integer k corresponding to the winding number of 
the Poly akov loops J27). 

For finite temperature and B = this model produces a crossover for both tran- 
sitions. Fig. [3] displays the condensates as functions of the temperature, and the 
critical temperature is defined by the change curvature in the curves. This occurs 
simultaneously for the chiral and deconfinement transitions within this model. 

At zero temperature, the Polyakov loop variable does not play a role. The pres- 
ence of a magnetic field enhances the chiral symmetry breaking, increasing the 
value of the chiral condensate, in line with the phenomenon of magnetic cataly- 
sis ll6l [8l l9l [181 17TI . This is shown in Fig. [4] It also deepens the minimum of the 
potential as B is increased, as illustrated in the same figure for several values of the 
magnetic field 

The dependence of the chiral condensate on the magnetic field is approximately 
linear, as shown in Fig. [4] This is in line with results from chiral perturbation theory 
|[T4l [T31 [TBI . Recent lattice results |47ll seem to deviate from a linear behavior for 
large B, growing faster, in better qualitative agreement with results from PNJL ll30l . 
However, for larger values of B, all model calculations seem to deviate from the 
lattice data, whereas for very small B they all quantitatively agree l47l . 
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Fig. 3 Expectation values of the order parameters for the chiral and deconfinement transitions as 
functions of the temperature. The filled circles represent the c-condensate, and the empty circles 
stand for the expectation value of the Polyakov loop. In this plot the condensates are dimensionless. 
Extracted from Ref. 1271 . 

Turning on the temperature, one can investigate the effects of the magnetic field 
on the thermodynamics and phase structure of strong interactions as captured by 
this model description. In the confining sector, the strong magnetic field affects the 
potential for the expectation value of the Polyakov loop via the intermediation of 
the quarks in three ways (27) : (i) the presence of the magnetic field intensifies the 
breaking of the global Z3 symmetry and makes the Polyakov loop real-valued, as 
shown in Fig. 12 (ii) the thermal contribution from quarks tends to destroy the con- 
finement phase by increasing the expectation value of the Polyakov loop; (iii) on 
the contrary, the vacuum quark contribution tends to restore the confining phase by 
lowering the expectation value of the Polyakov loop. 

In fact, the vacuum correction from quarks has a crucial impact on the phase 
structure. If one disregards the vacuum contribution from the quarks, as was done in 
Ref. [25), one finds that the confinement and chiral phase transition lines coincide. 
Moreover, in this case an increasing magnetic field lowers the equivalent chiral- 
confinement transition temperature. On the other hand, the inclusion of the vacuum 
contribution from quark loops in a magnetic field modifies completely the picture: 
confinement and chiral transition lines split, and both chiral and deconfining critical 
temperatures become increasing functions of the magnetic field. Both scenarios are 
shown in Fig. [5] which exhibit the full calculation of the phase diagram from the ef- 
fective potential within the PLSM (/ effective model. The vacuum contribution from 
the quarks affects drastically the chiral sector as well. Our calculations also show 
that the vacuum contribution seems to soften the order of the phase transition: the 
first-order phase transition - which would be realized in the absence of the vacuum 
contribution - becomes a smooth crossover in the system with vacuum quark loops 
included. 
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Fig. 4 Upper: the expectation value of the (dimensionless, a = E,v) condensate as a function of 
the magnetic field. Black dots are obtained from the PLSMq and the orange line is the linear fit. 
Lower: effective potential for the condensate at zero temperature for several values of the magnetic 
field B. Extracted from Ref. (27). 



The modifications produced by strong magnetic fields over strong interactions 
seem very exciting, bringing new possibilities for the phase diagram: affecting the 
nature of the transitions, splitting different coexistence lines, possibly exhibiting 
new phases, increasing the breaking of Z3, and so on. As discussed in the Introduc- 
tion, the second scenario has been also found in other effective models containing a 
chiral and a Polyakov loop sector [30 39], as well as in preliminary lattice simula- 
tions [44 1 . However, lattice simulations of magnetic QCD with physical masses and 
fine grids have shown that both critical temperatures actually go down for increasing 
B, saturating for very large fields ll46l . an unexpected behavior that is very different 
from the scenario depicted above. 
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Fig. 5 Phase diagram in the B-T plane. Upper: without vacuum corrections: the critical temper- 
atures of the deconfinement (the dash-dotted line) and chiral (the dashed line) transition coincide 
all the way, and decrease with B. Lower: with vacuum corrections: the critical temperatures of the 
deconfinement (the dash-dotted line) and chiral (the dashed line) transition coincide at B = and 
split at higher values of the magnetic field. A deconfined phase with broken chiral symmetry ap- 
pears. The vertical line is the magnitude of the magnetic field that expected to be realized at LHC 
heavy-ion collisions [J). Extracted from Ref. (27). 



This leads us to consider of a much simpler model that, yet, seems to contain 
the essential ingredients to describe the behavior of the deconfming line, and pro- 
duces results that are in qualitative agreement with the lattice: the magnetic MIT bag 
model [|40l . As will be clear in the discussion, the subtraction procedure in renor- 
malization is subtle (which can be seen as the choice of the renormalization scale) 
but can be guided by known physical phenomena and lattice results. 
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4 Magbag - the thermal MIT bag model in the presence of a 
magnetic background 

In the MIT bag model framework for the pressure of strong interactions, one needs 
the free quark pressure. As seen previously, the presence of a magnetic field in the 
z direction affects this computation by modifying the dispersion relation to 



I = 0,1,2,... being the Landau level index, s = ±1 the spin projection, / the fla- 
vor index, and qf the absolute value of the electric charge. Loop integrals are also 
affected as presented previously 12511511 . 

Since it has been shown that only very large magnetic fields do affect signifi- 
cantly the structure of the phase diagram for strong interactions ll24l l25l l29l l30l l44l 
[46), we can restrict the free quark pressure to the limit of very high magnetic fields, 
where it is possible to simplify some analytic expressions. 

It is crucial to realize, however, that the lowest Landau level (LLL) approxima- 
tion for the free gas pressure is not equivalent to the leading order of a large magnetic 
field expansion. For the zero-temperature, finite-B contribution to the pressure, the 
LLL is the energy level which less contributes in the limit of large B; the result be- 
ing dominated by high values of t. Nevertheless, the equivalence between the LLL 
approximation and the large B limit remains valid for the temperature-dependent 
part of the free pressure (as well as for the propagator), simplifying the numerical 
evaluation of thermal integrals l49l . 

The free magnetic contribution to the quark pressure has been considered in dif- 
ferent contexts (usually, in effective field theories [25 28 29 37 , 72]) and computed 
from the direct knowledge of the energy levels of the system, Eq. (1201 . The exact 
result, including all Landau levels, has to be computed from 



where the first term is a clearly divergent zero-point energy and the other one is the 
finite-temperature contribution for vanishing chemical potential. Since COm grows 
with B, the largest the I labeling the Landau level considered the larger the zero- 
point energy term becomes, being minimal for the LLL, corroborating the previous 
discussion. Thus, in the limit of large B, the LLL approximation is inadequate here. 
The decaying exponential dependence of the finite-temperature term on (Q( s f, on 
the other hand, guarantees that the LLL dominates indeed this result for intense 
magnetic fields. 

To obtain a good approximation for the large B limit of the free pressure, we 
choose to treat the full exact result and take the leading order of a x/ = mj-/2qfB 
expansion in the final renormalized expression. Let us then discuss the treatment 
of the divergent zero-point term. Despite being a zero-temperature contribution, the 
first term in Eq. (|2~TT > cannot be fully subtracted because it carries the modification 



G>lsf{k z ) = kj + mj + q f B(2£ + s + l) = k 2 +M| 



(20) 




(21) 
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to the pressure brought about by the magnetic dressing of the quarks. Using dimen- 
sional regularization and the zeta-function representation, which is also a type of 
regularization, for the sums over Landau levels and subtracting the pure vacuum 
term in (3 + 1 ) dimensions, one arrives at: 

2 

p l = SL(*/*)1c'(-M/)+^/-4)W+t 



2tt 
1 

12 



2/e + log(A 2 /2 ?/ fi) + i; 



(22) 



where a pole ~ (qfB) 2 [2/e] still remains. This infinite contribution that survives the 
vacuum subtraction can be interpreted as a pure magnetic pressure coming from the 
artificial scenario adopted, with a constant and uniform B field covering the whole 
universe (analogous to the case of a cosmological constant). In this vein, one may 
neglect all terms ~ (qfB) 2 and independent of masses and other couplings (as done, 
e.g. in Refs. If25ll28l[37l ). concentrating on the modification of the pressure of the 
quark matter under investigation. This can be seen as a choice for the renormaliza- 
tion scale after the renormalization of a ~ F^F^ term representing the magnetic 
field, as discussed, e.g. in Ref. l37l . We will come back to this point in the sequel. 

The final exact result for the free pressure of magnetically dressed quarks is there- 
fore 



{q f B) 7 
y 2k 2 



fj? _ y 



C (-l,x f ) - C (-1, 0) + - (x f - x 2 ) \nx f + ~-L 



2n 2 I ~ 



(23) 



In Refs. (25J1281137J, the constant £' (-1,0) = -0.165421... was not subtracted. In 
the case of pions, however, the full subtraction ensures that magnetic catalysis, i.e. 
an enhancement of chiral symmetry breaking, at zero temperature 0E|5]|T8][7T), is 
realized. On the other hand, if this term is left, the pion contribution to the effective 
potential for the chiral condensate at large magnetic fields will eventually raise the 
minimum instead of lowering it. 

In the limit of large magnetic field (i.e. Xf = m 2 /(2qfB) — > 0), we obtain 



Pq_l^By (q f B) 



2 r 



2n 2 



x/ln. 



f 2n 



dkAn 



(24) 



Adding the free piece of the gluonic contribution and the bag constant SS, the 
pressure of the QGP sector in the presence of an intense magnetic field reads: 



/&3P = 2(A£-1) 



7t 2 T 4 

90 



(25) 
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It is clear that, for V eB much larger than all other energy scales, the pressure 
in the QGP phase increases with the magnetic field, which seems to favor a steady 
drop in the critical temperature with increasing B that would lead to a crossing of 
the critical line with the T = axis at some critical value for the magnetic field. 
However, the behavior of T C (B) also depends on how the pions react to B, so that 
the outcome is not obvious. 

In the confined sector, which we describe by a free pion gas, one may follow 
analogous steps in order to compute the contribution from the charged pions, which 
couple to the magnetic field, arriving at 



An 2 



eB „ f 



(26) 



where M 2 e = m\ + (21+ \ )eB and x K = m z t /(2eB). In this final expression all 
terms ~ (q/B) 2 and independent of masses and other couplings were subtracted, 
as discussed before. Notice that the spin-zero nature of the pions guarantees that 
all charged pion modes in a magnetic field, differently from what happens with the 
quark modes, are B-dependent. So, in the large magnetic field limit the thermal in- 
tegral associated with iz + and n~ is exponentially suppressed by an effective mass 
> (rr^ + eB), as was also noticed in Ref. |25l . and can be dropped. In this limit, we 
have 

^+^- 1 ^ B -^C (M) (-l,l/2)x 7r , (27) 

where 1,1/2) = -ln(2)/2 = -0.346574- • •. Neutral pions do not couple 

to the magnetic field and contribute only with the usual thermal integral l52l . 

As before, for s/eB much larger than all other scales, the pion pressure rises with 
the magnetic field, as a consequence of the subtraction of all terms that are inde- 
pendent of temperature, masses and other couplings in the renormalization process, 
which renders the pressure positive. Differently from the quark pressure, however, 
the B = pion pressure takes over for temperatures of the order of the pion mass, 
which is not small and always enlarged by the presence of a magnetic field (given 
its scalar nature). Moreover, for large T, the magnetic pion pressures converge to 
(1/3) of the B = pressure, since n° is the only degree of freedom that contributes 
thermally for large B. 

Each equilibrium phase should maximize the pressure, so that the critical line in 
the phase diagram can be constructed by directly extracting T C (B) from the equality 
of pressures. It is instructive, nevertheless, to consider a plot of the crossing pres- 
sures, as shown in Figure|6] The figure shows, as expected, a decrease in the critical 
temperature (crossing points) as B is increased due to the corresponding positive 
shift of the QGP pressure. However, T c seems to be saturating at a constant value. 
One can see that the critical pressure (crossing point) goes down, but then it bends 
up again due to the increase in the pion pressure with B. This combination avoids a 
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Fig. 6 Crossing pion gas and QGP pressures as functions of the temperature for different values 
of the magnetic field: eB = (black, solid, right-most), 20m 2 ,, 40m 2 ,, 60m 2 K (magenta, dash-dotted) 
and eB = 100m^ (gray, solid, left-most), where in, = 138 MeV is the vacuum pion mass. Extracted 
from Ref. (40l. 



steady and rapid decrease of the critical temperature, as becomes clear in the phase 
diagram shown in Figure[7] In fact, inspection of the zero-temperature limit of Eqs. 
d24b and d2Tb shows that there is no value of magnetic field that allows for a vanish- 
ing critical temperature. 




eB /m 



Fig. 7 Phase diagram in the presence of a strong magnetic field. We also keep the T c (B = 0) point. 
The blue square represents a very conservative estimate for the maximum value of eB expected to 
be achieved in non-central collisions at the LHC with the formation of deconfined matter. The ar- 
row marks the critical temperature for eB r; 2\0m 2 K (73|, expected to be found at the early universe. 
Extracted from Ref. 1401. 
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The phase diagram in the plane T — eB shows that the critical temperature for 
deconfinement falls as we increase the magnetic field. However, instead of falling 
with a rate that will bring it to zero at a given critical value of eB, it falls less and less 
rapidly, tending to saturate at large values of B. Remarkably, this qualitative behav- 
ior agrees quite well with the most recent lattice results with physical masses l46lFi 
As discussed in the Introduction, previous models ll25l l27l [30l |3T1 l32l [391 . have 
predicted either an increase or an essentially flat behavior for the deconfinement 
critical line as B is increased to very large values. The same was true for previous 
lattice simulations P4l . which could be reproduced by the authors of Ref. l46l by 
increasing the quark masses to unphysical values. 

The renormalization procedure in the presence of a constant and uniform mag- 
netic field seems to be very subtle and crucial for the phenomenological outcome 
for the phase structure. B-dependent, mass-independent terms survive pure vacuum 
(B = 0) subtraction and have to be subtracted either in an ad hoc fashion |40l or by 
including a background field counterterm associated with a term ~ F^F^ repre- 
senting the magnetic field 071 . The latter brings a renormalization scale and, upon 
an appropriate choice, reproduces the former. Subtracting all purely magnetic terms 
in the pressures seems to be the appropriate choice since: (i) one guarantees that 
the pion pressure grows with increasing magnetic field at zero temperature, which 
is consistent with the well-known phenomenon of magnetic catalysis; (ii) lattice 
simulations usually measure derivatives of the pressure with respect to temperature 
and quark mass, and do not access derivatives with respect to B, so that purely in- 
dependent terms are not included in their results; and (iii) the effect of a purely mag- 
netic contribution to the pressure would only shift the effective potential as a whole. 
In particular, there would be no modification on relative positions and heights of 
different minima that represent different phases of matter. 

The qualitative success of the description of the deconfinement transition in the 
presence of an external magnetic field in terms of the MIT bag model suggests that 
confinement dynamics plays a central role in guiding the functional behavior of T c . 
In this case, a large N c investigation of the associated magnetic thermodynamics 
seems appropriate. 



5 Large N c 

Lattice QCD calculations [ 74 1 show that the deconfinement phase transition of pure 
glue SU(N C ) gauge theory becomes first order when N c > 3 1731176117711781 with a 
critical temperature given by l79l 

lim = 0.5949(17) + (28) 



Of course, our description necessarily predicts a first-order transition, as usual with the MIT bag 
model, and our numbers should be taken as rough estimates, as is always the case in effective 
models. 
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where a ~ (440 MeV) 2 is the string tension. The thermodynamic properties of pure 
glue do not seem to change appreciably when N c > 3 l80l[8l1 . which suggests that 
large N c arguments may indeed capture the main physical mechanism behind the 
deconfinement phase transition of QCD. 

It has been shown in Ref. |50l that the deconfinement critical temperature must 
decrease in the presence of an external magnetic field in the large N c limit of QCD, 
provided that quarks exhibit a paramagnetic behavior. Assuming that Nf/N c <C 1 
and m q = 0, the only contribution to the pressure of the confined phase that enters 



at 



is given by the vacuum (B = 0) gluon condensate c^Nf.a . The gluon and 
quark condensates change in the presence of a magnetic field lfT31[T6 82 1 but these 
modifications are negligible in the large N c limit. Besides, the gluon contribution to 
the deconfined pressure is blind to the magnetic field. 

On the other hand, the quark contribution is affected by the magnetic field and 
has the form 



P q uark{T,eB) ~ N c Npairs{Nf) T f quar k{T / y/o ,eB/T 



(29) 



with Npai rs {Nf)/Nc <C 1 being the number of pairs of quark flavors with electric 
charges { (N c — 1 ) /N c , — 1 /N c } in units of the fundamental charge. Only the largest 
(~ iVj) charge in each pair contributes to leading order in Nf/N c . Notice that the 
function f qua rk is positive definite and must increase monotonically with T for a 
fixed value of eB until it goes to 1 in the high temperature limit T s/a, eB. Thus, 
one should expect that the critical temperature as a function of the magnetic field, 
T c (eB), must decrease with respect to the pure glue critical temperature, T c , by 
an amount of 0(Nf/N c ). This can be seen directly by equating the pressures at T c , 
which yields 11501 



TcjeB) 1/4 



T c (eB)\ _ c 2 (N P airs,eB) 



CSB 



(30) 



where we defined 



C 2 {N pairs, eB) =c 



1 N pairs {Nf) C « SB 



fquark 



(()) 



eB 



Vc ' r c (0)2 



N c 



SBJg 



I if! 

Jg lue \ ^/ct 



(31) 



Since C2(N pa i rs ,eB) < cq, one finds that T c (eB) /Tj ' < 1 by an amount ~ Nf/N c 
||50ll . Assuming that quarks behave paramagnetically for all values of B, then 
C2(Npai rs ,eB) < ci (Nf), its equivalent in the case withB = and Nf > 0, and T c (eB) 
is also lower than the critical temperature in the presence of Nf/N c flavors of mass- 
less quarks at B = ll50l . 

In a free gas implementation of the deconfined phase f„i m = 1 and, for very 
strong magnetic fields, f quar k ~ eB/T^ ll40l . so that the magnetic suppression of the 
deconfinement critical temperature goes like eBN pa irs/ {N c o). This simple imple- 
mentation in the limits of low and high magnetic fields provides a scenario in which 
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Fig. 8 Cartoon of the T c x eB phase diagram in the large N c limit, using the approximation of free 
deconfined quarks and gluons. The numerical value 0.59 shown in the plot was extracted from Ref. 
l79l . Extracted from Ref. (50). 

the slope in T c (eB) decreases for large fields, as illustrated in Fig. [8] which hints 
for a saturation of T c as a function of eB, as observed on the lattice 0461 and in the 
magnetic MIT bag model [ 4-0 1 , but cannot be obtained in a model-independent way. 



6 Conclusions and perspectives 

The investigation of the effects brought about by the presence of a magnetic back- 
ground on the thermal chiral and deconfining transitions is in its infancy yet. Never- 
theless, the promise of the outcome of a rich phenomenology in mapping this new 
phase diagram of strong interactions is concrete. 

First model calculations have revealed the possibility of modifications in the na- 
ture of the QCD phase transitions, and also the appearance of a new phase of strong 
interactions in the case of a splitting of the critical (chiral and deconfining) lines. 
Even if recent, more physical lattice simulations have drastically modified the initial 
picture, they have also shown that the magnetic background has a very non-trivial 
influence on strong interactions. For instance, the behavior of quark condensates at 
finite temperature is non-monotonic l47l . rendering well-established vacuum phe- 
nomena such as magnetic catalysis more subtle at finite temperature. 

The functional behavior of the critical temperatures still has to be understood 
more deeply. Although no model foresaw the fact that both, chiral and deconfining 
temperatures, decrease then saturate at a nonzero value according to the lattice 11461 . 
a posteriori the magnetic MIT bag model was successful to describe this behavior 
for deconfinement qualitatively BOl and seems to capture some essential ingredi- 
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ents. A model-independent analysis in the large-A^ ( limit of QCD also points to this 
behavior |50l , which is reassuring from the theoretical standpoint. 

Another key ingredient in building an understanding of the physics of the quark- 
gluon plasma under these new conditions, which can be relevant for high-energy 
heavy-ion collision experiments, the primordial quark-hadron transition and mag- 
netars is the standard perturbative investigation of magnetic QCD. The calculation 
of the pressure in thermal QCD to two loops in the strong sector using the full QED 
propagator in the lowest Landau level approximation is subtle but possible, as done 
originally in Ref. |49l . 

The computation makes use of the full magnetic propagator that was obtained by 
Schwinger 1831 . but can be cast in a more convenient form in terms of a sum over 
Landau levels as derived in Ref. l84l (see also Refs. |49l[85ll ). In particular, it has 
been shown in Ref. J49| that the chiral limit for the exchange diagram seems to be 
trivial for very large magnetic fields. Concretely, it can be written diagrammatically 
in the following compact form |j49l : 




which realizes the intuitive expectation that the nontrivial dynamics in an extremely 
intense magnetic field should be one-dimensional. Since gluons do not couple di- 
rectly to the magnetic field, their dispersion relation maintains its three-dimensional 
character, which effectively results in a "massive" gluon in the dimensionally- 
reduced diagram. In the end the exchange contribution to the pressure is essentially 
an average over the effective gluon transverse mass m\ = k\+k\ of the exchange di- 
agram in (1 + l)-dimensions with the Gaussian weight (qfB/2n)exp[— m\/2qfB]. 
Since the trace in the reduced diagram is proportional to mi, the chiral limit seems 
trivial ll49ll . A detailed analysis of the dependence of the pressure on the mass and 
temperature and a semiclassical interpretation of this result will be reported soon 

ma. 

The nature of the phase diagram of strong interactions in the presence of a mag- 
netic background is still open. Recent lattice data, especially when compared to 
effective model predictions, seem to indicate that confinement dynamics plays an 
important role in the phase structure that emerges and should be incorporated in any 
effective description. Comparison between lattice data with very different quark 
masses |44, 46, 47) also show that the dependence of the critical temperatures on 
this parameter is non-trivial: T c increases at the percent level for large masses l44l 
whereas it decreases appreciably for physical masses l46l . This competition be- 
tween the effects from the magnetic field and quark masses on T c was also found 
in the large-A^- QCD analysis of Ref. l50l . A more systematic analysis of this phe- 
nomenon on the lattice would be very helpful for the building of effective models. 
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